Interfacial tension and nucleation in mixtures of colloids and long ideal polymer coils 
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Mixtures of ideal polymers with hard spheres whose diameters are smaller than the radius of 
gyration of the polymer, exhibit extensive immiscibility. The interfacial tension between demixed 
phases of these mixtures is estimated, as is the barrier to nucleation. The barrier is found to scale 
linearly with the radius of the polymer, causing it to become large for large polymers. Thus for 
large polymers nucleation is suppressed and phase separation proceeds via spinodal decomposition, 
as it does in polymer blends. 



I. INTRODUCTION 

In earlier work Q] immiscibility in mixtures of colloidal 
hard spheres and long ideal polymer molecules was stud- 
ied. Long means that the root-mean-square end-to-end 
separation of the polymer molecules, Re, is larger than 
the diameter a of the hard spheres. A mixture of spheres 
and long polymers was found to demix at comparable 
number densities of polymer molecules and spheres, both 
densities scale with R E as l/R E a, for R E > a. This scal- 
ing comes directly from the leading order term in the sec- 
ond virial coefficient for the sphere-polymer interaction. 
The interaction and hence the virial coefficient must be 
extensive in the number of monomers for large Re and 
hence must scale as R\ for our ideal polymers. The re- 
quirement that it has the dimensions of a volume then 
imposes the scaling R%a as a is the only other relevant 
length in the problem: the monomer size is assumed to 
be much less than a and so is irrelevant. Once the mix- 
ture has demixed we have two coexisting phases: one 
with a high density of colloidal particles and a low den- 
sity of polymer molecules, and one with a high density 
of polymer molecules and a low density of colloidal par- 
ticles. There is an interface between these two coexisting 
phases. Here we determine the scaling of the interfa- 
cial tension 7 of this interface, and use it to show that 
when a mixed sample of polymer and colloid is prepared 
and quenched into the two-phase region, the dynamics 
of the separation into two phases starts off with spinodal 
decomposition not nucleation. The fact that the phase 
separation starts off with spinodal decomposition makes 
mixtures of the long polymers like polymer blends but 
unlike simple mixtures, e.g., mixtures of oil and water. 
Thus, we can apply much of what we have learned of 
spinodal decomposition in systems like polymer blends, 
to mixtures of colloidal spheres and much larger polymer 
molecules. 

The colloid-polymer interactions in and bulk thermo- 
dynamics of, mixtures of colloids and large ideal poly- 
mers have both been studied, see Refs. @-§, but this 
is the first study of the interfacial tension and nucleation 
of these systems. The opposite limit to that of interest 
here, i.e., where the polymer molecules are smaller than 




the colloidal spheres, has been considered extensively, see 
for work on the bulk phase behaviour and 
L3f for work on interfaces. The following two 
sections deal with the interfacial tension, and with nu- 
cleation. Throughout, the objective will be to determine 
the scaling of the behaviour with the ratio of the size of 
the polymer to that of the sphere. Also note that here 
the polymers are always ideal, mixtures of polymers with 
strong excluded volume interactions and spheres, are very 
different fSBl. 



II. INTERFACIAL TENSION BETWEEN THE 
DEMIXED PHASES 

The interfacial tension between the demixed phases, 
one colloid-rich, the other polymer rich, can be estimated 
using just dimensional analysis. The tension 7 is an en- 
ergy per unit area. It is obtained by multiplying the 
free energy per unit volume, which is kT/{R 2 E a) Q|, by 
the width of the interface. This width will be of the or- 
der of the polymer size Re- Thus, 7 ~ kT/RE<J- Note 
that the energy scale has to be the thermal energy kT 
as there are no other relevant energy scales in the prob- 
lem. The mixture is athermal, there are no attractive 
interactions or soft repulsions to provide another energy 
scale. The free-energy density is then of order kT times 
the number density, which is of order l/R E a for both the 
polymer molecules and the colloidal spheres when they 
demix. This is just a simple scaling argument so we con- 
firm it by determining the scaling of 7 within a standard 
square-gradient or Cahn-Hilliard theory for the interface 

We apply this theory to the system in the semigrand 
ensemble of Ref. |j| where the characteristic thermody- 
namic potential is the semigrand potential to, which is a 
function of the number density of colloidal particles, pc, 
and the activity z of the polymer molecules. As we are 
specifying the activity not the density of the polymer our 
system is equivalent to a single component system whose 
thermodynamic state depends on the density and on the 
activity of the polymer z; In z acts as an inverse temper- 
ature in the sense that the larger it is the stronger is the 



effect of the attractions. Thus we can apply the standard 
square-gradient expression for the interfacial tension of a 
single component system, which is [Qpj,Q5j], p2| 
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where 
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is the excess grand potential at a point, lu^ and p^f* 
are the semigrand potential and density in either one of 
the coexisting phases. The superscript a = C,P for the 
colloid-rich and polymer-rich phases respectively, pc is 
the chemical potential of the colloid. The interface is 
normal to the x-axis. The coefficient K of the gradient 
term is assumed to be density independent. The equi- 
librium profile is obtained by minimising Eq. ([!]). Then 
standard manipulations enable a simpler expression for 
the equilibrium interfacial tension to be derived [ pl| 
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If we require that the functional Eq. (Q) be consistent 
with linear response theory p2]| we obtain an expression 
for the coefficient k of the gradient term 
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where C2(r; pc-, z) is the direct correlation function of the 
fluid of colloidal hard spheres in the presence of polymer. 
For our systems the most basic assumption is to use the 
low density approximation to the direct correlation func- 
tion. This replaces c with the Mayer f function for the 
effective sphere-sphere interaction in the presence of the 
polymer p2y . For two spheres with centres separated 
by less than a, the interaction energy is infinite and the 
Mayer f function equals —1. For separations larger than 
a the only interaction is that due to the polymer. This 
interaction is known || , and is long-range and weak thus 
we linearise the Mayer f function. Adding this altogether 
we obtain 
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We now return to Eq. (Q) for the interfacial tension 
and determine its scaling with Re- We note that the 
density difference p^ — ~ l/R E a, k scales as 

given by Eq. (0), ~ kT/R E a and the polymer activ- 
ity is of order the polymer number density z ~ l/R E a. 
Putting this all together we see that we recover the scal- 
ing 7 ~ kT / Re<J obtained earlier by dimensional anal- 
ysis. Also, from Eq. (Q) we see that the characteristic 

length-scale for the interface must be (n/^) 1 ^ 2 p^f 1 ~ Re, 
as we assumed earlier. Earlier work by Vrij |llf |, and by 
Brader and Evans jl2| on the interfacial tension between 
demixed colloid-rich and polymer-rich phases when the 
colloid and polymer were of comparable sizes, R E ~ a, 
found that, as expected. 7 ~ kTja 2 ~ kT/R 2 E . This is 
consistent with experimental findings p3| . 

The interfacial tension 7 will be of order kT/R E cr 
only if we are not too close to the critical point of 
the polymer-colloid demixing. In general we have 7 = 
(kT / Re<j)s(z / z c — 1), where s is a scaling function and 
z c is the polymer activity at the critical point. We have 
been assuming that we are not very close to the criti- 
cal point, i.e., that z/z c — 1 = O(l), and for these val- 
ues of its argument s = 0(1) and we return to 7 being 
of order kT / Re<J- However, as the critical point is ap- 
proached, z/z c — 1 <C 1, we have that the scaling func- 
tion s(x) = sqx^ for x -C 1, where sq is a dimensionless 
constant and p is the (positive) critical exponent of the 
interfacial tension . The interfacial tension tends to 
as the critical point of demixing is approached, and near 
the critical point it varies as a power law. See the review 
of Widom j24j] for an excellent introduction to interfaces 
near critical points. Sufficiently close to the critical point 
the scaling of the interfacial tension will be dominated by 
fluctuations and then the exponent p will take the value 
for the Ising model in three dimensions, p = 1.26 [Ts|,p^|. 
However, for very long polymers Re 3> c the effective in- 
teraction is long-ranged and long-range interactions sup- 
press fluctuations and make the system mean- field-like. 
The mean- field value of the exponent fi is 3/2 Which 
value of the exponent, Ising or mean-field, is observed is 
determined by whether or not the Ginzburg criterion is 
obeyed or not; see Ref. or any introduction to criti- 
cal phenomena for a definition of the Ginzburg criterion. 
Note that Eq. (||), belonging as it does to a mean- field 
theory, will yield an interfacial tension which tends to 
with an exponent p — 3/2, its mean-field value. 



the ideal polymer induces an attraction which decays as 
1/r for separations less than the radius of the polymer 
and roughly exponentially beyond this. Putting our ap- 
proximate c into Eq. (Q) we obtain an estimate of this 
coefficient 
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III. NUCLEATION AND OTHER 
FLUCTUATIONS 



Now consider a single phase mixture of spheres and 
polymer quenched into the two-phase coexistence re- 
gion. For definiteness assume that the single phase is 
the polymer-rich one. Then in order for the second, 
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colloid-rich, phase to form and coexist with the polymer- 
rich one, this second phase must form. The dynamics of 
the formation of a new phase fall into two broad cate- 
gories: nucleation then growth, and spinodal decomposi- 
tion. See Refs. [^8| for an introduction to the dynam- 
ics of first-order phase transitions. For example a mixture 
of simple liquids such as water and an alcohol phase sep- 
arate via nucleation of the new water-rich or alcohol-rich 
phase followed by growth of the nuclei, whereas polymer 
blends phase separate via spinodal decomposition. Here 
we show that for large ideal polymers and spheres, nu- 
cleation becomes very difficult, so mixtures of large ideal 
polymers and much smaller spheres will start to phase 
separate via spinodal decomposition. 

The rate of nucleation N n can be estimated using clas- 
sical nucleation theory, see the book of Debenedetti 
for a comprehensive discussion, see also Refs. 0,^3] . N n 
is the number of nuclei crossing the barrier per unit time 
per unit volume. The classical nucleation theory expres- 
sion for the rate N n is 



N„ 



Texp(-AF*/kT), 
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where T is an attempt frequency per unit volume, gen- 
erally slowly varying, and AF* is the free energy barrier 
that must be crossed in order for a new phase to nucle- 
ate. The variation in the rate is generally dominated by 
that in AF* so we focus on this. The free energy barrier 
comes from the free energy needed to form a microscopic 
droplet of the new phase, here a colloid-rich phase. This 
droplet is the nucleus of the new phase. Within classical 
nucleation theory the free energy of formation of a mi- 
croscopic droplet is the sum of two terms, a bulk term 
and a surface term, 



AF = -ttF 3 *„ + 47ri? 2 7 , 
where R is the radius of the droplet and 



* n = *>o>P) 



w(pc) 



/ («) 

PC(P C 



(9) 



Pc) (10) 



is the difference between the grand potential inside the 
nucleus and the grand potential of the phase in which the 
nucleus forms, pc is the density of colloid in the phase in 
which the nucleus forms, and p^f 1 is the density of col- 
loid inside the nucleus. So long as we do not approach 
the spinodal too closely, we can express as a Taylor 
expansion in chemical potential, around the chemical po- 
tential of the colloid at coexistence [i co . Truncating the 
Taylor series after the linear term, we get 



~ (fie - Pco) , 

as ty n = at coexistence, with 



\dpcj MC=Mc 
the derivative of ^ n at coexistence 
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The barrier is given by the free energy of the droplet 
whose free energy is highest, which occurs when the two 
terms in Eq. (g) are comparable, F 3 5' n ~ R 2 ^. Thus, 
R ~ 7/(*'(Mc - Pco)), and 
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From Eq. (|l2|), scales as \/R 2 E a. Using this scaling 
together with that of 7, we find that at the top of the 
barrier the size of the nucleus is of order Re, for jic — Pco 
not too much less than kT. Using these same scalings in 
Eq. ([l3]), we obtain the principle result of this section, 
the scaling of the nucleation barrier AF*, 
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the barrier scales as Re I a and so increases as the size 
of the polymer molecules relative to that of the colloidal 
spheres increases. The factor in parentheses is a dimcn- 
sionless measure of the supersaturation: how far we are 
into the two-phase region. 

For sufficiently large Re of the polymers the nucleation 
barrier will become so large that nucleation cannot occur. 
The mixture will be metastable up to very close to the 
spinodal pp|, p5|-|27[| . Thus, the mixture will only start 
to demix when quenched beyond the spinodal, where the 
phase separation will start with spinodal decomposition. 
This is precisely analogous to polymer blends and sys- 
tems of particles in which the particle-particle attractions 
are long ranged. In these systems the nucleation barrier 
scales as TV 1 / 2 and as r 3 , where TV is the length of the 
polymer and r is the range of the attraction [^tJ. The 
phase transition dynamics of systems of polymers and of 
particles with long-range attractions, were studied exten- 
sively in the early 1980s by Binder and Klein and their 
coworkers pE-27|. Many of the conclusions of that work 
also apply to the demixing of mixtures of hard spheres 
and much larger ideal polymers. 

Finally, we note that our finding that nucleation is sup- 
pressed is equivalent to saying that our mixture satis- 
fies the Ginzburg criterion for the irrelevance of fluctua- 
tions p9| | . Essentially, nucleation is a fluctuation so when 
fluctuations are weak nucleation is suppressed and vice 
versa, again this was found for polymer blends/particles 
with long-range attractions [^7[^|. For our mixtures, 
the Ginzburg criterion is essentially that the root-mean- 
square (rms) fluctuations in the number of colloidal 
spheres or of polymer molecules, in a volume R E , are 
much less than the mean number in that volume. The 
volume R E is the volume over which a pair of spheres or 
polymer interact (the sphere-sphere interaction is medi- 
ated by the polymer molecules, and the polymer-polymer 
interaction is mediated by the spheres). The rms fluctua- 
tions scale as the square root of the number of spheres in 
a volume R E , which is ~ (Fb/ct) 1 / 2 , whereas the mean 
number scales as Re jo. Thus for large values of the ratio 
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of the sizes, Re/<?, the rms fluctuations in the number 
of spheres (or polymer molecules) inside the interaction 
volume is a small fraction of the mean number. Fluctua- 
tions about the mean are small and so mean-field theory 
applies and nucleation, which is a fluctuation, has a very 
high free-energy cost. 
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